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We analyze the AC Stark shift of the Cs microwave atomic clock transition theoretically and 
experimentally. Theoretical and experimental data are in a good agreement with each other. Results 
indicate the absence of a magic wavelength at which there would be no differential shift of the clock 
states having zero projections of the total angular momentum. 
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I. INTRODUCTION 

The present definition of the unit of time, the second, 
is based on the frequency of the microwave transition 
between two hyperfine levels of the Cs atom. Recently, 
it has been realized that the accuracy and stability of 
atomic clocks can be substantially improved by trapping 
atoms in optical lattices operated at a certain "magic" 

wavelength pig. At this magic wavelength, both clock 

levels experience the same AC Stark shift; the clock 
frequency becomes essentially independent on trapping 
laser intensity. 

This effect was demonstrated 0, 0, [H] for optical clocks 
using strontium atoms. An extension of this idea to mi- 
crowave clocks with alkali-metal atoms Rb and Cs was 
considered in Ref. [6[ . A multitude of magic wavelengths 
for the hyperfine transition was identified. Unfortunately, 
detailed analysis presented below shows the conclusions 
of that paper to be erroneous: there is no magic wave- 
length for Cs, at least for clock levels with zero projec- 
tions of the total angular momentum Mp on the quan- 
tizing magnetic field. In a separate paper Q we analyze 
the case of circular light polarization and Mp ^ levels 
and demonstrate that the AC shift can be eliminated by 
an appropriate "magic angle" choice of the direction of 
the magnetic field with respect to the light propagation. 

This paper presents a detailed theoretical analysis of 
the frequency shift of a microwave clock involving a hy- 
perfine transition. The analysis requires calculation of 
the differential polarizability involving third-order ex- 
pressions, quadratic in the field strength and linear in the 
hyperfine interaction. Evaluation of the resulting expres- 
sions is carried out using relativistic many-body theory. 
The second part of the paper reports measurements of 
the clock shift at two laser wavelengths. The results of 



the calculations are in a good agreement with the exper- 
imental measurements. 



II. THEORY 

Here we follow the formalism of the quasi-energy 
states reviewed in the context of laser-atom interaction 
in Ref. 8]. We start with considering the AC Stark 
shift in the second order of perturbation theory 

(quadratic in the electric field) and then extend the for- 
malism to the higher-order AC Stark shift 5E^ 2+1 ^ which 
takes into account the hyperfine interaction (HFI). The 
latter shift appears in the third order of perturbation the- 
ory and is quadratic in the field amplitude and linear in 
the HFI. An important part of the analysis involves the 
tensorial expansion of the shifts in the scalar, vector, and 
tensor parts. 

We are interested in transitions between two hyper- 
fine components of the same electronic states. Below we 
employ the conventional labeling scheme for the atomic 
eigenstates, \n(IJ) FMp), where I is the nuclear spin, 
J is the electronic angular momentum, and F is the to- 
tal angular momentum, F = J + I. Mp is the projec- 
tion of F on the quantization axis and n encompasses 
the remaining quantum numbers. Since the clock transi- 
tions involve the same electronic state, we will also use a 
shorthand notation \F,Mp). For example, the Cs foun- 
tain clock involves transitions between hyperfine levels 
\F = 4, Mp = 0) and \F' = 3, M' F = 0) of the 6s 1/2 
ground electronic state. 

Under the influence of the laser each clock level is per- 
turbed. The clock frequency is modified by the difference 
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in the perturbed energies 
<Ji/ Stark (w L ) = 



(1) 



At the "magic frequency", this AC Stark clock shift 
would vanish. 



A. Second-order dynamic response 

This section introduces notation and reviews deriva- 
tion and tensorial analysis of the conventional second- 
order dynamic Stark shift. We demonstrate that for 
states of total electronic angular momentum J = 1/2 
and for Mp = levels (or for linear polarization) the 
second-order AC Stark shift of the clock transition van- 
ishes. 

Consider a traveling electromagnetic wave of an arbi- 
trary polarization, 

£ = ^£ L e e -*M- fcz ) + c.c. , 

with the complex polarization vector parameterized by 
an angle 9 



i = e x cos (9 + ie y s'm0. (2) 

The parametric angle 9 may be related to the degrees of 
linear, / = cos 29, and circular, A = sin 29, polarization. 
Notice that the quantizing axis z is chosen along the 
propagation vector k of the laser. The field amplitude 
£l is related to the intensity of the laser as II = j^£li 

or in practical units I L [sgj] « 1.33 x (£ L [^]) 2 . 

In the dipole approximation, the coupling can be rep- 
resented as (h.c. is a hermitian conjugate) 

V E i (t) = -£-T> = --£ L e ■ De' lu)t + h.c. 

Application of the Floquet formalism (dressed states) 
yields the second-order AC shift of the atomic energy 
level a 
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, E a -(E b -u) Y E a -(Eb + uy 



where v = —\£l £ ■ D , ip a and ip b being the station- 
ary atomic states with unperturbed energies E a and E b , 
respectively. 

Now we proceed to the conventional reduction of the 
polarizability into a sum over irreducible tensor opera- 
tors. Introducing the resolvent operator (Hq is the un- 
perturbed atomic Hamiltonian) 



R Ea (w)= (E a -H 



(2) 

we may recast the shift as an expectation value SEa 

(i£) 2 (lP a \6 m H|Va), With 



Oei (w) 



{e-T>) i R Ea (w)(e-D) 

+ {£-T>)R Ea (-u) (e-D)*. 



The order of coupling of the operators may be changed 
Oei M 



e [(-^y 

K =0,1,2 

x {i* <g> e} K ■ {D ( 
+ {e* <g> e} K ■ {D i 



R Ea (w) B} K 
R Ea {-w)T>} K ] 



leading to the conventional decomposition into the scalar 
(K = 0), vector (K = 1), and tensor (K = 2) terms. 

Here we employed {e (g) £*} KM = (-1) {£* ® e} KM 
and the fact that e and D are rank 1 tensors. The 
Mk component of the compound tensor of rank K com- 
posed from components of the tensors Ak x and Bk 2 (of 
rank K x and K 2 ) is defined as {A Kl (g> B K2 } KM — 

J2m 1 m 2 C KiM 1 k 2 m 2 A k 1 m 1 Bk 2 m 2 , where C Ki] £ K:iM2 
are the Clebsch-Gordan coefficients. The gener- 
alized scalar product is defined as [Ak ' Bk) = 

Ea/ k (-1) Mk Akm k Bk -Mk- 

Using the Wigner-Eckart theorem, a matrix element 
between two atomic states may be expressed as 

(FMf\Oei (w) \FM' f ) = (-1) A 'E(- 1 ) M 

K=0,l,2 n 

x {i' 



F K F 
-M F fi Mp 



F—Mf 



with the reduced polarizabilitics 



a. 



(K) 



( W ) 



(nF\\{D® R EnF {oj)D} k 

+ (-lf{D(g> R EnF {-u)D} K \\nF) 

VW}(-i) k+2F E{f FF> 

pi L 

x Y,( nF \\ D \\ n ' F ')( n ' F '\\ D \\ nF ) 

n' 

i 



+ (-i) 



£"„f — E n > p i +uj 

K 1 



E n l 



E 



n F 



(3) 



Here we have used the shorthand notation [A"] = (2K + 
1). The matrix element may be simplified further us- 
ing specific parametrization, Eq. @, of the polariza- 
tion vector. Explicitly,{e* <g> e} 00 = {e* 8) e} lM = 
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Finally, the AC Stark energy shift reads 



SE. 



(2) - - I -Sr 

nFM F — \ 2 



<f (w) + A a a nF (w) 



3M| - F (F + 1) 
2F(2F- 1) 



~2F 



(4) 



with the conventional scalar, vector, and tensor polariz- 
abilities 



<F («) 



1 1 



(°) c \ 

1 1 2F 
"TlyiF] 7F(F + 1) 



V6 



2F(2F- 1) 



(2F - 2) 



(5) 

1/2 



(2F + 3)! 



In general, there is an off-diagonal Mp = M' F ± 2 optical 
coupling involving the tensor part of the polarizability. In 
practice, a quantizing B-field is applied along the prop- 
agation of the laser wave, and as long as the off-diagonal 
coupling is much smaller than the Zceman intervals, it 
can be disregarded. 

Since the dipole matrix elements do not couple to the 
nuclear degrees of freedom, the dependence of the re- 
duced polarizabilities on I and F may be be factored 
out as 



45? h 



(-i) 



J+I+F+K 



J F I I _(K) , , 

F J k r a nJ 



where the quantities a^P (u>) are the reduced matrix el- 
ements in the \nJMj) basis, 

atfpiu) = (nJ\\{D® R EnF (w) D} K 

+ (-lf{D® R EnF {~u)D} K \\nJ). 

These quantities do not depend on either I or F. 

With this factorization, we can make important com- 
ments specific to the case of J = 1/2 (e.g., ground state 
of alkali-metal atoms such as Rb and Cs). Due to the 
angular selection rules, the tensor contribution (expec- 
tation value of the rank 2 tensor) vanishes and the only 
contributions come from the scalar and vector parts. As 
the vector part of the energy shift is proportional to Mp, 
for Mp = clock levels only the scalar contribution re- 
mains. The vector contribution also vanishes for the case 
of linear polarization (A = 0). 

The next important step is to demonstrate that the 
scalar shift does not depend on F. In other words, there 
is no clock shift at the second order. Indeed, for the scalar 
term, a s nF (u) = -^-^_ a W( w ) = ^_i_a<°>( W ). 

This result holds for an arbitrary J. 

This result has a very simple explanation. As it was 
pointed out above, for the case of linear polarization of 
laser light and J = 1/2 we have only a scalar contri- 
bution to the polarizability. This contribution does not 



depend on the orientation of the quantization axis (ex- 
ternal magnetic field). Let us consider the case when 
the quantization axis is directed along the laser electric 
field (along the linear polarization vector e). From the 
symmetry of the problem it is obvious that the electron 
states | J z = 1/2 > and \J Z = —1/2 > have exactly equal 
quadratic shifts. The hyperfine states |F, F z > are linear 
combinations of these electron states multiplied by nu- 
clear states. Since both components, \ J Z = 1/2 > and 
\J g = —1/2 >, of any hyperfine state have the same shift, 
all hyperfine states have the same shifts, i.e. the differ- 
ential polarizability is equal to zero. To have a non-zero 
differential polarizability one has to include the hyperfine 
interaction. Note that the inclusion of the magnetic po- 
larizability does not change this conclusion. If the lattice 
laser frequency is in optical range, the magnetic polariz- 
ability contribution is suppressed by an additional factor 
fi 2 B /D 2 ~ a 2 w (1/137) 2 and may be neglected. 

To summarize, we arrive at the conclusion that for 
J = 1/2, Mp = clock levels (or for linear polariza- 
tion) the second-order AC Stark shift is zero. Since the 
calculations of Ref. [f| were limited to this second-order, 
their conclusions are erroneous. This is also shown in the 
Appendix without using irreducible tensor algebra. 



B. Non-trivial effect of the hyperfine interaction 

In the previous section, the HFS interaction has served 
the role of an "observer" , as it only defined the coupling 
scheme. In particular, we find that the Mp = levels 
of the hyperfine manifold attached to the J = 1/2 levels 
are shifted identically - at that level of approximation 
any laser wavelength is "magic", i.e., the clock transition 
remains unperturbed by any laser field. The non-trivial 
effect arises when we take into account the dynamic (as 
opposed to the observer) role of the HFS interaction. 

Formally, this effect appears in the third-order double 
perturbation theory with two laser and one HFS inter- 
actions. We build the perturbation theory in terms of 
combined interaction 

V = V hls + V E x (*) ■ 

The convenience of the Floquet formalism is that we 
may immediately employ the conventional formula for 
the third-order energy correction 



VabVbcVcc 



-V aa 2^ 



(0) 



&#a (El 



(0) 
b 



eP 



Here a, b, c are dressed atomic states, i.e. \a) = \ip a )e muJt , 
with n representing the number of photons (n could be 
both negative and positive). The scalar product, in addi- 
tion to the conventional Hilbert space operational defini- 
tion includes an averaging over the period of oscillation of 
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the laser field. Explicitly, after the time averaging (now 
a, b, and c are atomic states) 



(a|tWRE. (0) vi?is a (w)v f |a) 

+ (a|F h fsi?B a (())«* JZjj. (-w)«|a), 

+(a|i; t i? £ ; a (-w) 14fe-RB a (-w)v\a) 

" (^fs) Qa (<o|«t (i? Ba H) 2 «|a) 
+ (aMi? Ba (-u;))V|a)' 



T„(w) 
S Q (w) 



Here T a (w), C a (w), and B a (oj) stand for top, center, and 
bottom position of the HFS interaction in the respective 
diagram. The term O a (cj) describes the normalization 
term. The relevant diagrams are shown in Fig. [TJ 



1/ 



(a) 



(b) 



(c) 



(ci) 



A 



(e) 



FIG. 1: Contributions to the dynamic polarizability a(cu). 
Diagram (a) represents the second-order contribution arising 
from two photon interactions (wavy lines). Diagrams (b-e) 
represent the additional effect of the hyperfme interaction 
(capped line) and correspond respectively to the third-order 
top, center, bottom, and normalization terms as described in 
the text. 
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The interaction of the electron with the nuclear mag- 
netic moment /i reads 



V hls 



(m-t«), 



where is the rank 1 irreducible tensor operator act- 
ing in the electronic coordinates with components 



(f) 



47T£ 



cr 



where a stands for the Dirac matrices and (f ) are 
the normalized vector spherical harmonics. In the formu- 
las below we require the reduced matrix element of the 
nuclear moment operator in the nuclear basis, (J||/*||J). 
It is related to the nuclear magnetic ^-factor as 



(7|| M ||/) = (2/+l)(2/ + 2) 



5M« 



H n being the nuclear magneton. 

Carrying out the angular decomposition similar to the 
second-order analysis of the preceding section we find 
that the expressions for the shift, Eq. (j4|), remain the 
same with the reduced polarizabilities o£w (ui) redefined 
as 



Q 



The third-order rank K = 0, 1, 2 corrections are given by 
the sum over contributions from diagrams of Fig. [TJ 

Pi K P M = 2$p (w; T) + ^ ( w; C) + (3<£) ( W ; O) . 
Explicitly, 



piV(^)=[F]mj:(-i) j+ ^{iy F }{(]jj }{f J F J F a jiis^), 



J a J b J, \ F F K 



J J Ji 

J a Jt 1 \cffi b {nJ,w), 
1 1 K 



(® o) = (-D 2J+1 [F] VW\ { F Yj } { J F J F * } E { H i } oSf «) 

Here we introduced the reduced sums 

Tjfj 6 (nJ,w) = (7|H|/)E E {nJ||T( 1 )||n a J a )(n a J a ||£)||n 6 J 6 ){n 6 J 6 || J D||nJ) 



l l 

E-E a E-E b +u; 



(6) 
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Cg 6 (nJ,u;) = (I\\n\\I) J2( n J\\ D \KJa)(n a Ja\\r {1) \\n b J b )(n b J b \\D\\nJ) 
0™(»J, M ) = (/ll M ||/)(nJ||T! I >||nJ)^{ n J||D|l»„J„){ n „J„|lDl|nJ) (— -J -j + (-1)" (<j - -u) ) .(8) 



Notice that the angular momenta of the intermediate 
states J a and J b are fixed. 



C. Numerical evaluation 

To perform the calculations we use an ab initio ap- 
proach which has been described in detail in Ref. . In 
this approach high accuracy is attained by including im- 
portant many-body and relativistic effects. 

Calculations start from the relativistic Hartree-Fock 
(RHF) method in the V N ~ 1 approximation. This means 
that the initial RHF procedure is done for a closed-shell 
atomic core with the valence electron removed. After 
that, the states of the external electron are calculated 
in the field of the frozen core. Correlations are included 
by means of the correlation potential method [13]. We 
use the all-order correlation potential X which includes 
two classes of the higher-order terms: screening of the 
Coulomb interaction and hole-particle interaction (see, 
e.g., [HI for details). 

To calculate X we need a complete set of single-electron 
orbitals. We use the B-spline technique [12] to construct 
the basis. The orbitals are built as linear combinations 
of 40 B-splines of order 9 in a cavity of radius 40as- The 
coefficients are chosen from the condition that the or- 
bitals are the eigenstates of the RHF Hamiltonian Hq of 
the closed-shell core. The X operator is calculated with 
the technique which combines solving equations for the 
Green functions (for the direct diagram) with the summa- 
tion over complete set of states (exchange diagram) [ll[ ■ 

The correlation potential X is then used to build a 
new set of single-electron states, the so-called Brueckner 
orbitals. This set is to be used in the summation in 
equations ©, ©, and ©. Here again we use the B- 
spline technique to build the basis. The procedure is very 
similar to the construction of the RHF B-spline basis. 
The only difference is that the new orbitals are now the 
eigenstates of the Hq + X Hamiltonian. 

Brueckner orbitals which correspond to the lowest va- 
lence states are good approximations to the real physical 
states. Their quality can be tested by comparing experi- 
mental and theoretical energies. Moreover, their quality 
can be further improved by rescaling the correlation po- 
tential X to fit the experimental energies exactly. We do 
this by replacing the Hq + X Hamiltonian with Hq + AX, 



I 

in which the rescaling parameter A is chosen for each par- 
tial wave to fit the energy of the first valence state. The 
values of A are A s = 1, X p = 0.97, and A<j = 0.95. Note 
that these values are very close to unity. This means that 
even without rescaling the accuracy is good and only a 
small adjustment of X is needed. Note also that since 
the rescaling procedure affects not only energies but also 
the wave functions, it usually leads to improved values 
of the matrix elements of external fields. In fact, this is 
a semi-empirical method to include omitted higher-order 
correlation corrections. 

Matrix elements of the HFS and electric dipole opera- 
tors are found by means of the time-dependent Hartree- 
Fock (TDHF) method [13, El- This method is equivalent 
to the well-known random-phase approximation (RPA). 
In the TDHF method, the single-electron wave functions 
are presented in the form ip = ipQ+6?p, where ipo is the un- 
perturbed wave function. It is an eigenstate of the RHF 
Hamiltonian Hq: (Hq — eo)ipo — 0. St/) is the correction 
due to an external field. It can be found by solving the 
TDHF equation 

(H - e )S^j = -SeiPo - Fifo - SV N -^ , (9) 

where Se is the correction to the energy due to the ex- 
ternal field (Se = for the electric dipole operator) , F is 
the operator of the external field (Vhfs or — D • £), and 
§yN-i j s correc tion to the self-consistent potential 
of the core due to the external field. 

The TDHF equations are solved self-consistently for 
all states in the core. Then the matrix elements between 
any (core or valence) states n and m are given by 

^nlF + SV"- 1 ^). (10) 

The best results are achieved when %p n and ip m are the 
Brueckner orbitals computed with rescaled correlation 
potential AX. 

We use equation (flU)) for all HFS and electric dipole 
matrix elements in evaluating the top, center, bottom, 
and normalization diagrams (Eqs. ©,©,©), except for 
the ground state HFS matrix element in the normaliza- 
tion diagram where we use experimental data. The re- 
sults are presented in section HVl 
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III. EXPERIMENT 



IV. RESULTS AND DISCUSSION 



We measure the frequency shift of the Cs clock tran- 
sition (|F = 3,M F = 0) to \F = 4, Mi? = 0)) in- 
duced by a far detuned laser beam. A Cs fountain clock 
[3| is used. At each cycle ~ 10 6 atoms are loaded in 
an optical molasses and cooled to below 2/iK. Moving 
molasses launches the atoms upwards with a speed of 
4.1 m/s where they pass twice through a microwave cav- 
ity thereby realizing Ramsey spectroscopy. The Zeeman 
degeneracy is lifted by a 1.6 mG magnetic field aligned 
along the fountains axis. The detuned laser beam is a 
traveling wave beam also aligned on the axis of the foun- 
tain. The light polarization is linear with respect to the 
light propagation. The beam waist is larger than the 
11 mm diameter opening in the microwave cavity. This 
assures that all atoms passing through the opening and 
being detected experience the light. The light intensity 
averaged over 1 cm 2 is measured by a commercial pow- 
ermeter (Newport 840-C) before entering the fountains 
vacuum chamber. One intensity measurement is taken 
before and one after each one-day run. The intensity 
drift between the two is of the order of 1%, however the 
error of the light intensity experienced by each atom is 
rather high. This is due to our ignorance of the exact 
intensity distribution, the exact atom distribution and 
intensity losses in the vacuum window as well as parasite 
reflections inside the vacuum chamber. We estimate the 
intensity error as 20%. 

The frequency shift is measured by alternating the 
fountain's configuration every 50 cycles. The first con- 
figuration is the standard clock operation. The second 
configuration is identical to the first plus the laser beam 
opened during the Ramsey period. This assures that the 
atom preparation and cooling is not disturbed by the 
laser and that the atomic cloud is identical in the two 
configurations. Hence, we can assume that all other clock 
shifts, in particular collisions, are identical for the two 
configurations. The absolute frequency is measured for 
each configuration against a highly stable local oscillator 
exhibiting no significant drift during several hundred cy- 
cles. The frequency shift induced by the light calculates 
as the simple difference 



The shift of the clock frequency is given by (cf. Eq. (Q])) 



where 8a {lol) is the differential polarizability. The 
conversion factor between differential polarizability in 
atomic units and the ratio of the shift to laser intensity 
in practical units is given by 



^Stark 



[Hz] 



j fmWl 



-4.68 x 10~ 5 x Sa (w) [a.u.; 



We start by considering DC polarizabilities. In the 
static regime (u>l = 0), our calculations give 5a — 
1.82 x 10 -2 a.u., which translates into the commonly used 
DC Stark coefficient k s = -2.26 x 10" 10 Hz/(V/m) 2 . 
Notice that this value includes only the scalar part of 
the polarizability. This is in agreement with the most 
accurate experimental result [15| of k$ = —2.271(4) x 
10- 10 Hz/(V/m) 2 . 
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cd 




(a.u.) 



8v — 2/ w ith light — ^without light- 



The frequency shift is measured for a number of laser 
intensities. Two sets of measurements are taken for light 
at wavelengths of 532 nm and 780 nm. The averages of 
each set weighted by the statistical frequency uncertainty 
give a light shift of (-3.51 ± 0.7) x 10" 4 Hz(W/cm 2 )- 1 
for 532 nm and (-2.27 ± 0.4) x 10~ 2 Hz(W/cm 2 )- 1 for 
780 nm. The statistical uncertainty is negligible before 
the uncertainty on the light intensity. 



FIG. 2: (Color online) Differential polarizability of the Cs 
clock transition (Mf = M' F = 0) in the B || k configuration 
as a function of the probe laser frequency. Two experimental 
measurements (at 780 nm and 532 nm) are compared with 
theoretical predictions (solid curve) . Refer to Fig. |3]for better 
graphical resolution of the 532 nm experimental point. 

For the AC case, our calculated differential polarizabil- 
ity 8a (lul) for the cesium clock transition is presented in 
Fig. [2] as a function of laser frequency. Both values are 
given in atomic units. The two peaks correspond to the 
6s — 6p and 6s — 7p resonances. The graph never crosses 
zero, which implies no magic frequency. Experimental 



7 



results for two laser wavelengths are also shown (also see 
Fig. [3]) . Calculated and experimental relative frequency 
shifts are compared in Table U and found to be in agree- 
ment with each other. 

0.05 



0.04 



TABLE I: Comparison of the theoretical and experimental 
AC frequency shifts for the clock transition in Cs 



0.03 



3 



■§> 0.02 



0.01 




Q I I I I I I I I I I I I I I I I I I l_ 

0.08 0.085 0.09 0.095 0.1 
a, (a.u.) 



FIG. 3: (Color online) Same as in Fig. [2] in the region of the 
experimental point at 532 nm. 



A, [rim] 


w, [a.u.] 


Svl/IL; 


[Hz/mW/cm 2 ] 






Theor. 


Expt. 


780 


0.0584 


-1.95 x 10~' A 


-2.27(40) x 10~' z 


532 


0.0856 


-3.73 x 10~ 4 


-3.51(70) x 10" 4 



To summarize, we presented a comprehensive analysis 
of the AC Stark shift of the Cs microwave atomic clock 
transition. Theoretical analysis based on the second and 
third order perturbation theory is accompanied by mea- 
surements. Calculations and measurements are in good 
agreement with each other and indicate the absence of 
a magic frequency at least for the Mp = clock levels 
with zero projections of the total angular momentum on 
the quantizing magnetic field. 



APPENDIX 

Considering the complexities of working with the an- 
gular algebra, here we analyze Eq. (2) of Ref. [6j for the 
2nd-order dynamic Stark shift. Starting from their equa- 
tion we again show that there is no AC Stark shift of 
the clock transition frequency. Eq. (2) of Ref. [f| for the 
polarizability contains a summation over M' , which we 
manipulate 



W F' 1 F \ _ ,_ iV f'-m+pW F 1 F' \f F' 1 F 
2^ \ M' p —M ) ~~ K ' Z^lM -p —M' J I M' p —M 



M' x ' M' 



'«wf_ F M k F M )(i P R oi){y F <-> 



A' 



The expression (|A. 1|) is obtained from using the summation rule 12.1(5) of Ref. [T(|; we further obtain expression (|A.2[) 
by noting that only Q — terms are non-zero in the summation. 

Now we take the .F'-dependent part of (|A.2[) with the .F'-dependent part of Eq. (2) of Ref. [6j and take the summation 
over F' 
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The expression (|A.3[) is obtained from using the summation rule 9.8(6) of Ref. [if 
Combining the above results gives 
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Not surprisingly, the six-j symbols here are identical to 
the ones appearing in the previously derived polarizabili- 

ties a£p(cv). Hence, this makes the connection to scalar 
(K = 0), vector (K = 1), and tensor (K = 2) parts. 

First we focus on the case p — 0; this corresponds to 
linear polarization in the B || i geometry. For J = 1/2 
atomic states the tensor part (K = 2) is necessarily zero 
due to selection rules in the six-j symbols (this is the 
case regardless of polarization). Furthermore the vector 
part (K = 1) is zero due to the fact that the top row of 
the second three- j symbol sums to an odd number (or 
see (|A.4[) below, with p = 0). This leaves us with only 
the scalar part (K = 0) to analyze. In this case, the 
r.h.s. simply reduces to 1/(3 [J]) = 1/6. Thus we can 
conclude that for linear polarization of this geometry, 
the 2nd-order dynamic Stark shift is F-independent for 
,/ = 1/2 atomic states. 

For the B || k geometry, the linear polarization is re- 
garded as an equal mixture of a + (p = +1) and a~ 
(p = —1) circularly polarized light. Again the tensor 
part is necessarily zero for J = 1/2. For the vector part, 
we note the three- j symbol relation 



Thus, when we take equal mixtures of <j + and a~ light, 
the vector contribution drops out. Again, we are left with 
only the scalar part. Not surprisingly, we again obtain 
the result 1/(3[J]) = 1/6 when taking equal mixtures of 
<j + and er~ light. 

The above results then generalize to any geometry for 
linearly polarized light. 
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